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Stochastic resonance with applied and induced fields: the case of 

volt age- gated ion channels 



M. ShakedQ, Z. Schuss^ 



Abstract 

q^". We consider a charged Brownian particle in an asymmetric bistable electrostatic po- 

tential biased by an externally applied or induced time periodic electric field. While 
the amplitude of the applied field is independent of frequency, that of the one induced 
by a magnetic field is. Borrowing from protein channel terminology, we define the open 
probability as the relative time the Brownian particle spends on a prescribed side of the 
potential barrier. We show that while there is no peak in the open probability as the 
frequency of the applied field and the bias (depolarization) of the potential are varied, 
there is a narrow range of low frequencies of the induced field and a narrow range of 
the low bias of the potential where the open probability peaks. This manifestation of 
stochastic resonance is consistent with experimental results on the voltage gated Iks 
PQ . and KCNQ1 potassium channels of biological membranes and on cardiac myocytes. 
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cr 1 Introduction 



Our recent experimental findings show unusual non-thermal biological effects of a periodic 



\Q \ electromagnetic field (EMF) of frequency 16 Hz and amplitude 16 nT (nano Tesla) on the 

potassium current in human Jk s and KCNQ1 channels [Tj. More specifically, we expressed 
O ' the -^Ks channel in Xenopus oocytes and varied the membrane depolarization between -100 

^ . mV and +100 mV and measured the membrane potassium current. The current with applied 

EMF peaked above that without applied EMF at membrane depolarizations between mV 
and 8 mV to a maximum of about 9% (see Figures [1] and [2]) . A similar measurement of the 
potassium current in the KCNQ1 channel protein, expressed in an oocyte, gave a maximal 
increase of 16% at the same applied EMF and at membrane depolarizations between -10 
mV and -3 mV (see Figure [3]). Similar experiments with L-type calcium channels showed no 



response to the electromagnetic field at any frequency between 0.05 and 50 Hz. 

In a related experiment [2], we applied electromagnetic fields at frequencies 15 Hz, 15.5 
Hz, 16 Hz, 16.5 Hz and amplitudes of the magnetic field from below 16 pT and up to 160 nT, 
to neonatal rat cardiac myocytes in cell culture. In the range 16 pT - 16 nT, we observed 
that both stimulated and spontaneous activity of the myocytes changed at frequency 16 Hz: 
the height and duration of cytosolic calcium transients began decreasing significantly about 
2 minutes after the magnetic field was applied and kept decreasing for about 30 minutes 
until it stabilized at about 30% of its initial value and its width decreased to approximately 
Vo. About 10 minutes following cessation of the magnetic field the myocyte (spontaneous) 
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Figure 1: Ik s current in Xenopus oocytes with applied magnetic field of 16 Hz and 16 nT (red) 
and without (blue) 



Iks relative change - ELF EMF: B = 16nT, f = 16Hz 




Figure 2: The quotient of Jks expressed in Xenopus oocytes with applied magnetic field of 16 Hz 
and 16 nT and without (red to blue in Figure [T]) 
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Figure 3: The quotient of KCNQ1 expressed in Xenopus oocytes with applied magnetic field of 16 
Hz and 16 nT and without 

activity recovered with increased amplitude, duration, and rate of contraction. Outside 
this range of frequencies and magnetic fields no change in the transients was observed (see 
Figure EJ). When the stereospecific inhibitor of KCNQ1 and Jks channels chromanol 293B was 
applied, the phenomenon disappeared, which indicates that the Jk s and KCNQ1 potassium 
channels in the cardiac myocyte are the targets of the electromagnetic field, in agreement 
with the former experiment. The effect of changing the outward potassium current in a 
cardiac myocyte is to change both the height and duration of calcium transients, action 
potential, sodium current, as indicated by the Luo-Rudy model [3]. 

The specific response at 16 Hz may indicate some form of resonance or stochastic res- 
onance of a gating mechanism of open voltage-gated potassium channels (e.g., a secondary 
structure or mechanism) with time-periodic induced electric field. Since the induced electric 
field is too low to interact with any component of the Iks channel, we conjecture that the 
induced field may interact with locally stable (metastable) configurations of ions inside the 
selectivity filter We propose an underlying scenario for this type of interaction based on 
the collective motion of three ions in the channel, as represented in the molecular dynamics 
simulation of [I]. The configurations of three potassium ions in the KcsA channel is repre- 
sented in [J] in reduced reaction coordinates on a three-dimensional free energy landscape. 
In our simplified model, we represent the collective motion of the three ions in the channel 
as diffusion of a higher- dimensional Brownian particle in configuration space. An imitation 
hypothetical energy landscape with a reaction path (indicated in red) is shown in Figures [5] 
and [6j Projection onto a reaction path reduces this representation to Brownian motion on 
one-dimensional landscape of potential barriers (see Figure [7]). The stable states represent 
instantaneous crystallization of the ions into a metastable configuration, in which no current 
flows through the channel, that is, they represent closed states of the channel. There is also 
a pathway in the multidimensional energy landscape that corresponds to a steady 
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Figure 4: Cardiac cells, 4 days in culture, were exposed to magnetic fields of magnitude 160 pT 
and frequency 16 Hz for 30 min. Characteristic traces of spontaneous cytosolic calcium activity 
(A,B,C,D) and of electrically stimulated (1 Hz) cytosolic calcium activity (E,F,G,H). Times are 
measure in seconds from the moment of application of the magnetic field. 

current flowing in the channel, e.g., an unobstructed trough in the energy landscape. Tran- 
sitions from the latter into the former represent gating events. In our scenario the motion 
between closed states is simplified to one-dimensional Brownian motion, e.g., in a trough 
obstructed with barriers, while the interruptions in the current correspond to exits from the 
unobstructed trough into the obstructed one. Activated transitions over barriers separating 
two closed states in the obstructed trough (see Figure |Sj) affect the probability of transition 
from closed to open states. Stochastic resonance between two closed states may change the 
transition rates between them, thus affecting the open (or closed) probability of the channel 
(see Section @J. 

We investigate the stochastic resonance (SR) in our mathematical model of a Brownian 
particle in an asymmetric bistable potential with an induced electric field. The difference 
between this problem and that of the extensively studied SR with an applied periodic electric 
field [5], [6] is that according to Faraday's law (or Maxwell's equations), the amplitude of 
the induced field is proportional to the frequency of the applied magnetic field. While the 
traditional manifestation of SR is a peak in the power spectral density of the trajectory of 
the resonating particle, we consider its manifestation in the probability to be in one of the 
two meta-stable states. This measure of SR is ineffective for a symmetric potential, because 
this probability is 1/2 in the symmetric case and is independent of the applied periodic field. 
It is effective, however, in asymmetric potentials, for example, when a constant bias field 
depolarizes the membrane, as is the case in the above mentioned experiments. Note that 
in the second experiment the depolarization of the myocyte membrane is due to the action 
potential in the cell. In contrast, asymmetry of the potential can weaken SR with an applied 
field, as shown in [TJ, [8]. 

Our main results concern SR with applied and with induced external periodic forces. In 
the former case, which we view as a benchmark for our method of analysis, we find that 
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Figure 5: Hypothetical energy landscape of two ions in the selectivity filter. The reaction path is 
marked red. The straight segment in the trough may represent the open state in the channel 




Figure 6: Another view of the hypothetical energy landscape of two ions in the selectivity filter. 
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Figure 7: Profile of one-dimensional electrostatic potential landscape biased by a constant electric 
field 



2 parabolas profile 




Figure 8: A simplified version (see eq.Q) of the wells in Figured The wells at x\ and £3 
and the barrier at x<i are now at X\ = —2,x^ = 0.7 and Xi = 0. The constant bias in ([2]) is 
c = 0. 
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Normalized Open Probability 




Figure 9: P (u, c)/P o (0, c) with induced force A, A = 0.007, e = 0.029, x L = -2.4, = 
1.385 for < uo < 6, -0.1 < c < 0.1. Evidently, there is no SR. 




Figure 10: P (u, c)/P o (0, c) with induced force Aw, A = 0.007, £ = 0.029, x L = -2.4, a; R = 
1.385 

there is no SR as frequency and depolarization are varied, in agreement with known results 
[5] (see Figure [9]). In contrast, the probability to be on one side of the barrier in the case of 
an induced field peaks at a nearly fixed frequency in a finite window of depolarizations (see 
Figure [10]). We refer to this peak as stochastic resonance, though it may not be the usual SR 
phenomenon. The folding of the surface in Figure [9] into that in Figure [H)] seems to be due 
to the decrease in the amplitude of the induced field at low frequencies. This observation is 
consistent with the above mentioned experiments and seems to be new. 

To connect the above SR with the cardiac myocyte experiment, we use the Luo-Rudy 
model [3] of a ventricular cardiac myocyte of a Guinea pig. We express the manifestation 
of the above SR in the Ho dgkin- Huxley equations [9] as a change in the conductance of 
the Iks channel in the specific range of depolarizations at the resonant frequency of 16 Hz. 
We note that the Ik s is one of the delayed rectifier K + channels that are present in cardiac 
myocytes [TO] , in neuron cells [11], [12], and more, that is, it stays open long enough for its 
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Figure 11: Cytosolic potassium concentration [mM] vs time [msec] without SR (blue) and with 
SR (red) in the Luo-Rudy model. SPECIFY UNITS 



Intracellular calcium concentration - red : SR 
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Figure 12: Cytosolic calcium concentration [mM] vs time [msec] without SR (blue) and with SR 
(red) in the Luo-Rudy model. 



(secondary) gating to partially synchronize with the induced field. The SR-increased efflux 
of potassium (see Figure [TT]) shortens the action potential, and consequently lowers the peak 
of the cytosolic calcium concentration (see Figure [T2D, at the expense of increased sodium 
concentration (see [T3"j) . The shortening of the action potential leads to the shortening of 
the QT interval (see Figures HH [15]) [10] and was actually observed experimentally [13], 
P3]. These predictions of the SR modified Luo-Rudy equations are also new. In addition, 
we obtain from the SR modified Luo-Rudy model an increased conductance during the 
plateau of the action potential in the cardiac myocyte. This in turn shortens both the action 
potential and the cytosolic calcium concentration spike durations, lowers their amplitudes, 
increases cytosolic sodium, and lowers cytosolic potassium concentrations. These theoretical 
predictions are supported by experimental measurements. Specifically, these effects were 
communicated in [13], [H], as well as in our own measurements [2]. 



8 



ntracel lular sodium concentration - red : SR 



.1315 1.132 1.1325 1.133 1.1335 1.134 1.1345 



Figure 13: Cytosolic sodium concentration concentration [mM] vs time [msec] without SR (blue) 
and with SR (red) in the Luo-Rudy model. 




Figure 14: Action potential [mV] vs time [msec] without SR (Blue) and with SR (Green) in the 
Luo-Rudy model 
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Figure 15: Action potential duration [msec] vs time [msec] without SR (Blue) and with SR (Red) 
in the Luo-Rudy model. 

2 The mathematical model 



We consider the dimensionless overdamped dynamics 

d(f)(x, t) 



x 



dx 



4> x (x,t) 



(i) 



in the bistable time-periodic potential 

(f){x,t) = (c - A App i jInd ) smut)x + (f> (x), (2) 

where A^pi^d is the amplitude of the applied (induced) electric field and (j>o(x) is a fixed 
parabolic double well potential that consists of the two parabolas 

(x-x L ) 2 



4>o(x) 



[X - X R 



— 1 for x < 



1 for x > 0, 



(3) 



x 



R 



where x^ < < xr. The amplitude of the electric field induced by the time-periodic magnetic 
field B cos cut (B = const) is Aj n d = Au, where A = CB and C is the proportionality constant 
in Faraday's law. The linear term cx represents the membrane depolarization. This model 
can be considered the limit of the parabolic double well potential that consists of the three 
parabolas 
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for x > xg R , 



(4) 



10 




0.5 - 



-0.5 - 




, 1 1 1 1 [ 1 1 1 

5 10 15 20 25 30 35 40 



Figure 16: The deterministic trajectories (in dimensionless units) with uj — 1, A = 0.5, c = 
—0.05, xl = —2.4, xr = 1.385 are attracted to the periodic x^it) (blue lower curve), XM{t) 
(green middle curve) and xji(t) (red upper curve). 



where xl < —%s L < < xs R < xr and a > 0. The three parabolas connect smoothly 

at — x$ T and xs„, which implies the relationships xs T = %L . , x$ R = Xr , , 

1 + ax z L /2 1 + ax R /2 

lim axs L = —2/xl and lim axs R = 2/xr. The potential <ft(x,t) (see Figure ED has two 

a — >oo a — >oo 

periodic attractors, Xi(t) and XR(t) (see Figure [16]) and the separatrixo XM(t) = 0. The 
attractors are the stable periodic solutions of (CD), given by 

Xi{r) = at - AiCos(r + tpA, i = L,R (5) 



•; ^-Appl,Ind x 1 , 2 / Cxf \ 

Aj = — , (fij = arctan — ^— , a,- = x, . 6 

v/4 + jy' ^ v 2 ' 



When small white noise \/2ew(t) is added to the dynamics (JTJ, it becomes the stochastic 
equation 



± = -<^(x,t) + V2£w(t). (7) 

The trajectories of (J7|) spend relatively long periods of time near the attractors xi,(t) and 
xn(t), crossing at random times. The first passage time from xi(t) to xn(t) is defined 
as 

r L (t ) = inf{t > : x{t ) = x L (t ), x{t + t)= x R (t + t)} (8) 
and the mean first passage time is defined as 

f L = fj Q ®T L {t )dt , (9) 



3 In the model (Tj| of three parabolas the separatrix is XM^t) — ~ — ( g a^,a j [wcoswi + asinwi]. 
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where E denotes ensemble averaging over trajectories of ((Tj) and the period is 

T = — . 

The first passage time tr(£o) and the mean first passage time tr are defined in an anal- 
ogous manner. The fraction of time the random trajectory x(t) spends in the basin of 
attraction of x R (t), that is, the fraction of time that x(t) > XM{t), is the right probability 
P R (c, u, A Appl j nd , e), given by 

-nT poc 



pni poo 

P R {c,u,A App i Jnd ,e) = lim — / / p(x,t)dxdt 

pT poo 

= lim — / / p(x, t + nT) dx dt 

7WOO T J Q J S .„( t+nT ) 



1 

T 



Jx M {t+nT) 
T poo 

I p OQ (x,t) dxdt, (10) 

Jx M (t) 



where p(x,t) is the transition probability density function (pdf) of the random process x(t), 
generated by the stochastic dynamics (J7J) and Poo(x,t) = lim p(x, t + nT) is the periodic 

n— >oo 

pdf. We obtain in a similar manner 

pT px M (t) 

P L (c,uj,A App ij nd ,e) = — / p 00 (x,t)dxdt, (11) 

1 JO J-oo 

For small e, 

P L {c, u, A ApplJnd , e) « — — — . (12) 

tr + t l 

3 The Fokker-Planck equation 

The T-periodic pdf Poo(x, t) is the T-periodic solution of the Fokker-Planck equation 

dp{x,t) d 2 p(x,t) d[<j> x (x,t)p(x,t)] n/ + / 

— = e — — 1 — 1 for — oo < x < oo, < t < oo. (13) 

at ox 2 ox 

We construct a WKB approximation to Poo(x,t) for small e, 



exp | 


tf;(x,t,e)\ 
. e J 




/>oo 

/ exp • 

</ — oo 


/ */>(x,t,e) 
I e . 


| dx 



Poo(x,t) -5 J , (14) 



where ij;(x,t,e) is a T-periodic regular function of e. Expanding 

ip(x, t, e) = ip(x, t, 0) + eipi(x, t) + . . . , (15) 
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we find from large deviations theory that ip(x, t, 0) is the+ minimum of the integral 



I{x{-)){x,t) 



[x(s) + 4> x (x(s), s)} 2 ds 



(16) 



over all continuous trajectories x(-) such that x(0) = x. Setting r = ut, we write the 
Hamilton- Jacobi (eikonal) equation for the minimal values of I(x, r) in the domains x > 
and x < as 



^(x,r,0) = -(^) 2 (x,r,0)-- 

00 UJ 



c - A App i >Ind smT + 



2(x - Xi) 



.17 



<(x,r,0), (17) 



for i = L,R. The solution can be constructed in the quadratic form [T7] 

nx,r,o) = [x ~ Xi } T)]2 + ai (is) 

and the constants a, are determined from the Freidlin-Wentzell extremum principle [16J. 
According to this principle the local minima ip L (x,r, 0) and ip R (x,T,0) are joined into a 
global minimum function tp(x, r, 0) by the requirement that the steady state probability 
current across the separatrix % vanishes [TTj . 



/■27T 

J(0,r) = / [J L (0,T) + J R (0,r)}dr = 0, 
Jo 



(19) 



where the probability flux density is 

6^(0,7") 



./i(0,T) 



+ s (O,ry(O,r) 



6^(0, r) 



ax ax 
Using the WKB approximation (TH|) for jj(x, r), we find that the minimum condition is 

1 ftfir 



(20) 



4(0, r) = - 
Using (J2U in (HSJ, we find that 

Ja L -a R )/e _ 



-2x,(r) 



exp 



+ a, 



.17 



(21) 



4 J 


/" x L (r)exp< 








4 J 


£fl(r) exp < 


r 4^1 


> dr 


l J 



(22) 



It should be noted that the phases (pi and (p R may be disregarded in the integrals of equation 
fl22|) . therefore we set them to zero. Expanding the integrals in (122]) by the Laplace method 
for small e about the maxima of the integrands, at 



xi = x R (0) = a R - A R , x 2 = x L (ir) = a L - A L , 



(23) 
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we get 



e (a L -a R )/e 



where, according to (jHJ) and (pi = 




2irex 2 L „ . . ( x 2 l (tc) 



~%V>( ^ L {n exp I 2 



fl i fi exp 




a; 2 exp < \ 

-1A L x 2 I ex 2 L 



2 

I 3/ 1 



x l\hr? a^exp 



(24) 



3.1 The left probability P L ( 

To calculate the probability Pl(c, uj, AAppijnd, e), we use the WKB approximation (Til]) in 
( 1TT1) and evaluate the integrals by the Laplace method, as in (ED, to get 



-Pl(c, u;, AApplJnd, £J 



1 



1 _|_ Xjl c (ai,-aRVe' 



(25) 



-XL 



Using the result from (j23J) in ([25]) . we find that 

Pl(c, LU, AApplJnd, Z) = 



1 + 



X 



R 



V4 + xX 



^l| 3 V a/4 + x%u 2 V »i 



-x 2 



exp 



,'L i 



2 2 
&0C SOC 



(26) 



We normalize the frequency- dependent left probability Pl(c, u, AA pp ij n d, e) by the left prob- 
ability of the unforced dynamics P°(c, a;, AAppi,ind — 0, e). To calculate P°(c, o>, AAppijnd — 
0, e), we set A = in (1251) and obtain 



P°(c,u;,y4A pp z,/„d = 0,e) 



1 + 

X L 



2 - cx L 



2 - cxr 



exp 



c{x R - x L ) / c(x fl + x L 



4 



Note that (1271) is not other than (1121). where 



27T£ 



exp 



(27) 



(2f 



The MFPT in ( J28l) is the Kramers escape rate of a Brownian particle over a high sharp 
barrier xm [191. 
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4 Coarse-grained Markov model of secondary gating 



We consider the movement of a Brownian particle over two unequal barriers of heights 
A02i = 0(2:2) — <P( x i) an d ^043 — 0(2:4) — 0(2:3), respectively, such that A0 2 i <C A0 43 (see 
Figure Cj). Both x\ and £3 are closed states of the channel whereas £5 represents the open 
state (see Figure [5]). Our goal is to elucidate the influence of SR between the periodic force 
and the activation over the local small barrier A</> 21 at £ 2 , within the closed state, on the 
open probability of the channel. Specifically, when SR increases the time spent in the well 
at £ 3 relative to that at x±, the attempt frequency to cross the barrier at £ 4 into the open 
state £5 increases, thus increasing the open probability of the channel. More specifically, we 
evaluate the influence of SR on the mean closed time, that is, on the mean time spent in the 
wells at x\ and £3 prior to passage into £5 (which we denote ff 3 ). For that purpose, we can 
assume that £5 is an absorbing boundary. 

First, we note that steady state considerations can be applied in describing SR in the 
wells at Xi and x 3 . Indeed, we assume that 



which means that a transition over the barrier at £ 4 between the open and closed states occurs 
at a much lower rate than those over the barrier at £ 2 , between the two closed substates. 
In particular, the first inequality in (129]) means that there will be many transitions over 
the barrier at £ 2 before a transition occurs from £3 to £5 over the barrier at £4. Thus we 
confine our attention to transitions over the former and consider £5 to be an absorbing state, 
as mentioned above. The assumption of high barriers (the last inequality in ( 1291) means 
that a quasi steady state is reached in each of the wells before a transition over £ 4 occurs. 
Therefore the pdf of the quasi-steady state in each well can be represented by the principal 
eigenfunction and eigenvalue in that well, with absorbing boundary conditions. 

We coarse-grain the trajectory of the diffusion process x(t) into that of a continuous-time 
three state Markov jump process x(t), that jumps between x\ and £3 and is absorbed in £ 5 , 



The three state continuous-time Markov chain x(t) is not stationary due to the passage to 
the open state £ 5 . There are two time scales of passages: a short scale corresponding to the 
transitions between x\ and £3, and a long one for the transitions between £ 3 and £5. Due to 
the long time scale, the dynamics between the closed states [x\ and £3)is quasi-stationary. 
We assume it as a stationary dynamics in our analysis. 

The jump of the Markov process from x\ to x 3 occurs when x(t) reaches £3 for the first 
time after it was at £1, and so on. The Chapman-Kolmogorov equation for the transition 
probability matrix P t of the Markov process is [15] 



0(£ 4 ) - 0(X 5 ), 0(£ 4 ) - 0(£3) > 0(£2) - 0(^l) > 0(^2) - 0(^3) > B, 



(29) 




(30) 



Pt 



P t R, 



(31) 



where 



P t {i,j) = Pr{£(t) = j\x = i}, for 1 < i,j < 3 



(32) 



15 



and the elements of the instantaneous jump rate matrix R are 

-ri 3 ria 

R r 3 i -(r 3 i+r 35 ) r 35 



The stationary distribution 7r of the process is 

7T = [ 0, 0, 1 ] , 

because 25 is an absorbing state. 

Next, we calculate the time-dependent probability distribution 

p t = [ Pr{x(t) = xi}, Pr{x(t)=x 3 }, Pr{x(t) = x 3 } ] T . 

According to ( 13T1) and ( 1321) . p t the sum of elements in each column of the matrix P t 



Pr{x(t) = Xj} = Pr{x(t) = j I Xq = i} for 1 < j < 3 



and therefor satisfies the Chapman-Kolmogorov equation 



given by 



P t = R Pt 



Pt = e Po> 

where we assume that p is the stationary distribution of the chain x\ ^ x 3 , namely, 

r 3 i ri 3 



Po 



n 3 + r 31 na + r 31 



, 



We further express the vector p t as a linear combination of the eigenvectors {v\,V2,v 3 
the matrix R T , corresponding to the eigenvalues Ai, A2, A3, 



p t = a 1 e Xlt v 1 + a 2 e X2t v 2 + a 3 e Xat v 3 , 



where 



Ai 
A 2 



-^13^35 

s 



-s 



ri 3 r 35 
S 



A3 = 0. 
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and S = r 13 + r 31 + r 35 . Using the fact that lim p t = iz, we get that a 3 = 1. 



t— >oo 



The MFPTs r 3 and Ti are related to the exit rates r 31 and r 13 over a non-sharp boundary 
according to 



1 



T"3 



Tl 



2r 31 
1 



2r 



(42) 



13 



In order to find mean closed time rf 3 prior to the first arrival to 25, it is enough to 
consider the matrix 



-7*13 ^31 

ris -(ri 3 + r 35 ) 



because this time is determined by the first two elements of the vector p t , 

p 1 / = [ Pi{x(t)=x 1 }, Pr{x(t) = x 3 } f. 
An eigenvector expansion p]' 3 , similar to that in (138]) . is 

p]> 3 = a^V! + a 2 e X2t v 2 , 



(43) 



(44) 



(45) 



where Ai and A2 are given in ( |39l) and ( )40l) . respectively, and i>i and v 2 are the eigenvectors 



~ T 

of the matrix 



v\ = [r 3 iS,r 13 S - r 13 r 35 ] 



v 2 = [r 31 S,r 13 S - S 2 + r 13 r 35 ] T . 



Using the initial condition 



Pi" 



r 3 i 



ri3 



. r 13 + r 31 r 13 + r 31 J ' 



we obtain 



1 



a 1 



S 2 + r 13 r 35 

S 3 (r 13 + r 35 ) (r 13 + r 31 ) 2 ' 

1,3/ 



a: 2 



^13^35 



S 3 (r 13 + r 31 )' 



Setting P(t) = p]' 3 (l, 1) + Pt'^i^, 1) and substituting the values of a>i, a 2 , Vi and v 2 into 
151). we obtain 



P{t) = e Xlt - e X2t 



(ri 3 r 35 ) 2 
(ri 3 + r 3 i) 4 ' 



Hence 



r°° 1 

rt 3 = E[rl 3 ]= / P(t)rft~ — 
Jo l A i| 



S 



^13^35 ^35 V r 13 



(46) 



(47) 
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Using (j4"2"|) in fH7|) and setting P 3 fi = — , we obtain that 

T"3 + T\ 

^HK^- (48) 

We further coarse-grain the trajectories of the process into that of a telegraph process 
with two states, closed state (corresponding to x\ and 23) and open state (corresponding to 

c — > o. 

Denoting r D as the mean first passage time from x 5 , using equation (148]) we get 

-Popen = c = J J - • (49) 

T o "T T l,3 ^ 4. 

° r 35 P 3 R 



Applying the theory proposed in 13.11 to the dynamics between the closed states x\ and x 3 , 
the SR effect increases P R , by using a negative depolarization, c . According to fH9|) an 
increase of P R causes to an increase of P ope n- 



4.1 High barrier approximation to r% 5 , r$i, 
We consider the autonomous stochastic differential equation 

dx = -(f)'(x)dt + V2edw (50) 

x (0) = x. 

The transition rates between the wells are the probability fluxes in the direction of the tran- 
sition at the top of the barrier. Thus, denoting by 3>i(x) and A, the principal eigenfunction 
and eigenvalue in well i (i — 1, 3), we have 

n 3 = -8^X2), r 3i = £$3(22), r 35 = -£$ 3 (x 4 ). (51) 

To calculate the fluxes, we have to construct the eigenf unctions &i(x), which are the 
solutions of 

e®"(x) + [<j)'(x)$i(x)]' = -\i$i(x) for - 00 < x < x 2 (52) 

$i(x 2 ) = 0, $i(x)^0 for (53) 

£$3(2;) + [(j)' (x)$ 3 (x)]' = -A 3 <l>3(a;) for - x 2 < x < x 4 (54) 

$ 3 (x 2 ) = 0, $3(24) = 0. (55) 
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The asymptotic structure of the eigenfunctions is given in [19] as 

/~2~ PUJ 2 (x-X 2 )/y^£ 



2 ru 2 (x-x 2 )/^e 

$i(z) ~ -K l e~ m/£ \J- J e~ z /2 dz for x < x 2 (56) 

/2 [" ru 2 {x-x 2 )/Ve 
- / e~ z /2 dz - 1 for x 2 < x < x 4 , (57) 

where uOi = ^\<p"(xi)\, (i = 1, 2, 3, 4) and 

/ x 2 \/2tt r Xi \/2tt 

Q^dx*- — e-^ Xl)/e , Af 3 = <S> 3 (x)dx~- — e -* (xa)/e . (58) 
-oo U l Jx 2 ^3 

According to (JUD and (|56|)-(l58l), 



ri3 ^ glgg e -[0(xaW(xi)]/ e rgi _ ^3f^ e -[0(x a )-^(x 3 )]/e ^ _ ^4 e -[0(x4)-fl» 3 )]/ 
7T 7T 7T 

which are Kramers' rates for the corresponding barriers [101. 



5 Effect of SR in the Luo-Rudy model of cardiac my- 
ocytes 

The Luo-Rudy model [3] describes ionic concentrations and cardiac ventricular action poten- 
tial by a system of 21 ordinary differential equations. It reflects the guinea-pig electrophysi- 
ology by detailed Hodgkin-Huxley models of ionic currents. The most significant currents are 
the slow Jk s and rapid Ik t delayed rectifier potassium currents, a time-independent potas- 
sium current, a plateau potassium current (ultra-rapid Ik™-)) a transient outward current, 
fast and background sodium currents, L- and T-type calcium currents, a background cal- 
cium current, calcium pumps, sodium-potassium pumps, and sodium-calcium exchangers. 
In addition, the model describes Ca 2+ handling processes, that is, calcium dynamic release 
from the sarcoplasmic-reticulum and from the calcium buffers troponin, calmodulin, and 
calsequestrin. 

The stochastic resonance described above changes the open probability of the Jk s channel, 
and therefore it affects its conductance. To incorporate this effect into the Luo-Rudy model, 
we modify the Hodgkin-Huxley equation for the Jk s cur rent- volt age relation by shifting the 
stationary open probability of the channel in the Luo-Rudy model [3], 

p ° {y) = rVi. 5 ^ (60) 

1 + exp 

to 



16.7 



Po(V) = ? V + 4 . 12 ^ . <«) 



1 + exp 



16.7 
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Figure 17: Resonant increase (red) of 20% in the open probability of the Iks channel, and normal 
regime (blue) 

which imitates the experimentally observed shift (see figured]). This changes the channel 
conductance GksPo^Y) i n the Luo-Rudy model (Gk s is the open channel conductance) to 
<5ks-Po(^)' which changes, in turn, the membrane potassium current (Iks), averaged over 
many channels, to [3] 

(Iks) G Ks P*(V)(V - E Ks ) as t ^ oo (62) 

(Eks is the reversal potential of the channel). The effect of this modification of the Luo- 
Rudy model is shown in Figure [UJ The duration of the action potential is reduced and 
accordingly, the peak of the cytosolic calcium concentration is lowered (see Figures fl~2l) . as in 
the experiment described in the Introduction. On the other hand, sodium concentration is 
increased (see Figure [T3|) . The shortening of the action potential duration in the ventricular 
cardiac myocytes affects the QT interval in the electrocardiogram, which consists of a sum of 
several different action potentials created in the myocardium [TU] (see Figures [T ^ [To]) . These 
theoretical predictions are supported by experimental measurements. Specifically, these 
effects in vivo were communicated in [13], [TJ], as well as in our own in vitro measurements 

6 Conclusion and Discussion 

This paper tries to explain the results of the experiment of exposing human potassium Iks 
channels and cardiac myocytes, which contain these channels, to weak and slow electromag- 
netic fields. We offer a scenario of a new kind of stochastic resonance between the induced 
periodic field and the thermally activated transitions between locally stable configurations 
of the mobile ions in the selectivity filter. 

More specifically, since the induced electric field is too weak to interact with any compo- 
nent of the Iks channel protein, our model cannot describe the primary gating mechanism 



Open Probability of Iks channel 
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of a voltage gated channel. We therefore resort to a mathematical model, which postulates 
interaction of the induced field with configurations of the mobile ions inside the selectivity 
filter. These configurations may be much more susceptible to the weak induced field than 
any components of the surrounding protein, because the potential barriers separating the 
metastable configurations of the mobile ions can be of any height. 

According to our scenario, the observed resonance is due to the dependence of the induced 
electric field amplitude on frequency, in contrast to an applied external electric field with 
fixed frequency, which is known not to exhibit stochastic resonance with changing frequency 
and depolarization. In our theory the observed SR between two closed (or inactivated) states 
affects the open probability of the channel. 

Our model describes the dynamics of a Brownian particle in an asymmetric bistable 
potential forced by a periodic induced electric field. The analysis of this model is based on 
the construction of an asymptotic solution to the time-periodic Fokker-Planck equation in 
the WKB form. We evaluate the dependence of the steady state probability to be on one 
side of the potential barrier on the frequency, amplitude, depolarization, and noise intensity. 

Our main results are shown in Figure [10l which indicates that there is a peak in the 
open probability in a relatively narrow range of depolarizations and frequencies. We refer 
to this peak as stochastic resonance, though it is not be the usual SR phenomenon. This 
observation is consistent with the results of the Jk s channel experiment mentioned in the 
Introduction. 

Another result is the incorporation of the SR result into the Luo-Rudy model of cardiac 
myocytes. We found that the increased conductance of the Jk s channel reduces the duration 
of the action potential, the peak height of the cytosolic calcium concentration, in good 
agreement with the experimental results. The shortening of the action potential duration in 
the ventricular cardiac myocytes affects the QT interval in the electrocardiogram. 
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